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In one dimension the area law for the entanglement entropy is violated maximally by the 
ground state of a strong inhomogeneous spin chain, the so called concentric singlet phase 
(CSP), that looks like a rainbow connecting the two halves of the chain. In this paper we 
show that, in the weak inhomogeneity limit, the rainbow state is a thermofield double of 
a conformal field theory with a temperature proportional to the inhomogeneity parameter. 
This result suggests some relation of the CSP with black holes. Finally, we propose an 
extension of the model to higher dimensions. 
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I. INTRODUCTION 


Entanglement has become a very useful tool to study the structure of complex quantum states, 


such as the ground states (GS) of interacting systems [1]. Geometry and quantum structure are 
linked via the so-called area laws [2, 3|, which state that the entanglement entropy of a block B , 
calculated in the ground state of a local Hamiltonian, scales with the area of the boundary of the 
block. Area laws have been proved only in a few cases, such as ID gapped systems [4], and are 
violated in several interesting cases, such as conformal systems in ID [5H7|, where the entropy 
grows logarithmically with the block size, with a prefactor proportional to the central charge of 
the conformal field theory (CFT). These exceptions of the area law have received an increasing 
attention in recent studies [2, |f|. 

A strong violation of the area law takes place in an inhomogeneous free fermion model in ID 
where the hopping amplitudes between consecutive sites decay exponentially outwards from the 
center of the chain [lo|]. By tuning the exponential factor, the GS of the system evolves smoothly 
from a logarithmic law towards a volume law for the entanglement entropy between the left and the 
right halves of the chain [11], In the strong inhomogeneity regime, when the exponential factor is 
high enough, the GS is the product state of Bell-pairs symmetrically distributed around the center 
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of the system, as shown in Fig. |T] This valence bond state was termed concentric singlet phase 
(CSP) 10] or simply, rainbow state [hi]. The volume law for strong inhomogeneous chains can be 
easily understood from the rainbow picture. The entanglement entropy of half of the chain is given 
essentially by the number of bonds connecting the left and the right halves. However, for weak 
inhomogeneities the rainbow picture does not hold because the GS is a resonating valence bond 
state that satisfies a volume law plus logarithmic corrections fll]. Another type of inhomogeneity 
has been used to provide smooth boundary conditions which improve the convergence of the bulk 
properties of the GS [12|- Moreover, an exponential increase of the couplings has been used in 
a Kondo-like problem I bill and a hyperbolic increase for the study of the scaling properties of 
non-deformed systems [14, 15]. 

The aim of this paper is to understand the free fermion model in the weak inhomogeneity 
regime and its relation to the uniform limit given by a conformal field theory, namely a massless 
Dirac fermion with open boundary conditions. One might think that some scaling limit of the 
model would correspond to a perturbation of the underlying CFT. However, the perturbation 
cannot correspond to local operators added to the action, since they do not give rise to volume law 
entropies. Quite surprisingly, the solution of this puzzle is still provided by CFT: the GS is a sort 
of thermal state that satisfies a volume law, with a temperature related to the exponential factor 
of the hopping amplitudes. From this perspective, the appearance of a volume law in the GS of 
the model is not surprising at all since, after all, it corresponds to a thermal state. Yet, it comes 
as a surprise that the state remains pure. The understanding of this apparent contradiction will 
bring us to unexpected territories that we shall start to explore. 

The organization of the paper is as follows. After a general reminder of our model in section m 
containing an analysis of the strong inhomogeneity limit, we establish in section IIIII a continuum 
approximation in the vicinity of the homogeneous point, given by a deformation of the critical 
Hamiltonian. This deformation affects the single-particle wavefunctions and the low-energy ex¬ 
citations. In section El we study the entanglement entropy within the CSP. The continuum 
approximation of the previous section allows us to give an expression for the von Neumann and 
Renyi entropies of the left half of the system throughout the transition. Moreover, we show how 
the Renyi block entropies fit to the conformal expressions with varying coefficients. In section [V] 
we address the question: can the maximal area-law violation of the deformed hopping Hamiltonian 
be extended to higher dimensions? Indeed, we find that a natural extension of the concentric 
singlet phase can be found in 2D systems. The article ends with a summary of our conclusions 


and ideas for further work. Finally, in appendix[A] qubistic images m are shown to provide useful 
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information about the entanglement structure. 


II. MODEL AND NOTATION 


A. Strong Disorder Renormalization of the Hopping Model 


Let us consider a strongly inhomogeneous XX-model 

H X x = K i( a i a i+i + a ^ aV 


i u i+lh 


(1) 


i=l 


where the Ki > 0 are very different in value. We can apply the strong disorder renormalization 
group (SDRG) of Ma and Dasgupta [17] in order to obtain the GS. The renormalization prescription 
is to pick up the strongest coupling, K max , and to establish a singlet bond on top of it. Then, using 
second order perturbation theory, one finds the effective coupling between the two neighbours of 
the singlet, 

K L K R 


K = 


K n 


( 2 ) 


where Kr and Kr are, respectively, the left and right couplings to the maximal one. The renor¬ 
malization continues by choosing the next largest coupling and so on. Effective couplings may 
therefore emerge at long distances. The success of the SDRG scheme depends on the maximal 
coupling K max being always much larger than its neighboring values Kr and I\r. 

Using the Jordan-Wigner transformation, we can apply the same tools to study a fermionic hop- 

[ii, y, 


ping model in ID with the same features 

Jf. 


reading the couplings Ki as hopping amplitudes, 


L—l 


H = Ma+1 + h - c - > 


( 3 ) 


i= 1 


where c] creates a spinless fermion on site i, and is related to af via the non-local Jordan-Wigner 
transformation: 


c] = a f 


i —1 

n 

3 = 1 


er„- 


( 4 ) 


This non-local transformation points at a modification of the SDRG prescription to take into 
account the fermionic nature of the particles. Effective hoppings between non-contiguous sites 
are equal to the corresponding coupling in the XX model multiplied by a sign, J = (—1 ) np K, 
where uf is the number of fermions between the two sites [if] ] . This rule can be implemented with 
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a simple modification of the renormalization group (RG) prescription. Since a single fermion is 
always added at each RG step, 

j = -^4s. ( 5 ) 

Jmax 


This implies that the hoppings can be either positive or negative. When they are positive, a 
singlet-type bond is established between both sites, of the form |T~) oc 101) — 110). If the hopping 
is negative, the corresponding triplet-type anti-bond is established: |'F + } oc 101) + 110). Both types 
of bonds share many properties, such as the entanglement. They both represent different flavors of 
a Bell pair. The SDRG candidate for the GS of the system is a tensor product of bond or anti-bond 
states on the corresponding sites. Indeed, it can be written as a Fermi state: 

L/2 

i^) = n4io), (6) 

fc=i 


where d\. creates either a bond or an anti-bond on a pair of sites, i.e.: d\ oc c\ ± m. 


It was proved in 


181 ]. that the Hamiltonian Q has the following properties, for any values of 


the { Ji\. It presents particle-hole symmetry, which implies that for every single-particle eigenstate 
with energy e there is another eigenstate with energy — e, which is related by swapping the sign 
of the components of all odd sites. Thus, the ground state takes place at half-filling. Moreover, 
the occupation number of every site is equal: (n*) = 1/2. This is a non-intuitive result, given the 
inhomogeneity of Hamiltonian |1]), and it does not hold for excited states. 


B. The Rainbow State 


Let us describe the family of local Hamiltonians whose ground state approaches asymptotically 
the concentric singlet phase, also known as rainbow state, and give some heuristic arguments to 
explain the volume law for the entanglement entropy. 

Let us consider a chain with 2 L sites, which we will label using half-odd integers, n = 
±^, =t|, • • •, ± . On each site n, let c n and denote the annihilation and creation 

operators of a spinless fermion. The hopping Hamiltonian is given by [20] 


H = —— 


Ci c 


i_ — 

2 


T 3 

E 

__i 


Jn 


Cn^n+l T C_ n C_(j l _ 


+ h.c. 


( 7 ) 


where J n are the hopping amplitudes parametrized as (see Fig. [T]for an illustration) 
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| J 0 (a) = 1, 

I J„(a) = a 2n , n = ±,..., L - §. 


( 8 ) 


For a = 1, this is the uniform ID spinless fermion model with open boundary conditions (OBC). 
The model with 0 < a < 1 was introduced by Vitagliano and coworkers to illustrate a violation of 


the area law for local Hamiltonians [lcj] and was studied in much more detail in reference me 


Q 



Figure 1. Rainbow state showing the (—k, +k) bonds above the central link. Thus, the entanglement entropy 
between the left and the right halves of the chain is Llog2. 


In the a. —>• 0 + limit, the couplings become strongly inhomogeneous and, as argued in [ll], we 
can employ the SDRG described in the previous section. The largest hopping is the central one, 
Jo = 1, and gets renormalized to Jo = —a 2 , larger (in absolute value) than a 3 , which comes next. 
The values of J n are engineered so that this situation repeats itself for all steps of SDRG, so the 
bonds are established in a concentric way around the center, joining sites +k and —k, and giving 
rise to the aforementioned CSP, as in Fig. [I] Moreover, notice that the signs alternate. 

It is worth to notice the striking similarity between our system and the Kondo chain 21]. 


Indeed, let us divide our inhomogeneous chain into three parts: central link, left sites and right 
sites. The left and right sites correspond, in our analogy, to the spin up and down chains used in 
Wilson’s chain representation of the Kondo problem. In both cases, they form a system of free 
fermions, with exponentially decaying couplings. In the Kondo chain, notwithstanding, the central 
link becomes a magnetic impurity, which renders the full system non-gaussian. 


C. Strong Inhomogeneity Limit 

The a —>• 0 + limit of ([?]) , leading to the rainbow state, is singular: the Hamiltonian decouples 
in that limit, and only the central link survives. Let us consider a very small but non-zero a, and 
study the GS to first order in perturbation theory. The Hamiltonian is always free, so the GS is a 
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Slater determinant (0. The orbital operators, d\, can be expanded in terms of the local creation 
operators: 


4 = I>M 


kj 


(9) 


where are the wavefunction components for the single-body associated problem, i.e., eigenvectors 
of the hopping matrix. Let us propose the following form for them: 


_5 _3 _1 ,1 , 3 

2 2 2 ' 2 '2 


+ 5 
' 2 


"i/d = 

•0 2 = 

-i/i 3 = 


a 


a 1 1 a 

ala —a —1 —a 
1 a 0 0 a 1 


( 10 ) 


a 


Notice the sign alternation, due to the negative sign in the renormalization prescription Q. It 
is straightforward to check that all those ip k are eigenstates of the hopping matrix to first order in 
a. We can now define two families of states: the bonding and the anti-bonding creating operators, 
defined as: 

( 6 «)' = (4 + 4 )’ 

('’«)' = 71G - c ') ■ (11) 

0 + , the GS of Hamiltonian (jTj) can be written as the concentric singlet state 


In the limit a 
or rainbow state: 


\ Rl) - { b - L +h L -h) |0) - 


( 12 ) 


where sl = (—1) L . 


III. WEAK INHOMOGENEITY LIMIT: CONTINUUM APPROXIMATION 

The study of the weak inhomogeneity limit, a —» 1 _ , motivates the derivation of a continuum 
approximation of the Hamiltonian ()7|). This is obtained by expanding the local operator c n into 
the slow modes, iPr(x) and 4>l(x) around the Fermi points ±kp 

cz <? kFX Mx) + e- ikFX ^ R (x), (13) 

y/a 

located at the position x = an E (—£,£), where a is the lattice spacing and £ = aL. In the 
continuum limit, a —» 0 and L —» oo, with £ kept constant. At half-filling, kp = 7t/(2 a) is the 
Fermi momentum. 
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Equation (1131) is the familiar expansion used in the uniform case, a = 1, that will enable us to 
derive the numerical results found in 11]. Plugging Eq. (1131) into Eq. ([7]) one obtains 

rC 


tt m r A 

H ~ — / axe a 
2 J-c 


ipRdxi>R - {dxipR^R - ipLd x ilj L + {d x ipl)ip L 


/,t 


t, 


d' 


(14) 


where 


a = e~ h '\ 


(15) 


To derive equation GUI. we have assumed that the fields i/ Ri l(x) vary slowly with x. so that 
cross terms like (—1 ) x / a ^' R {x)^)L{x) can be dropped. We have also made a gradient expansion 
i/(x + a) ~ VK 3 ') + a d x ' l P(x), keeping only terms up to the first derivative. The Hamiltonian (1141) 
describes the low energy excitations of the original lattice Hamiltonian at half-filling. It is worth 
to mention that (1141) is a Hermitian operator, that is, = H, which is of course a consequence 
of the hermiticity of (|7|) . In the continuum limit we shall take h 0, with h/a kept constant, so 
that a —>• l~. 

The boundary conditions (BC) satisfied by the fields iPl,r(x) at x = ±C, can be derived from 
equation (fT3j) setting c ± ^ L+ i^ = 0 and taking a continuum limit that yields 

Vir(±£) = =rVt(±£)- (16) 


Then, integrating by parts, one can write Hamiltonian (1141) as 


h\x\ 


H ~ ia J dxe « 


V R d x ^R - i’ld x ^L ~ ^■sign(x)(V’] l V ; i? - V’LV’l) 


T. 


h 


/, t 


/,t, 


(17) 


The Fermi velocity is then given by vf = a, that we set equal to one by convention (similarly, we 
replace h/a —>• h). The single-body spectrum of the uniform model, that is, h = 0, can be easily 
found 


Em — 


7r(m + 1/2) 
2 L 


, m = 0 , ± 1 ,... 


(18) 


For the non-uniform model we have the Eqs. 


le 


—h\x\ 


d x T ^ sign(x) 


i/r,l{x) = ±Ei/) RtL (x), 


(19) 


whose solution is 


^r,l(x) = A RtL e h \ x ^ 2 


exp 


=F^sign(x) (e h \ x \ - l) 


( 20 ) 
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Notice that, in the limit h —>• 0, one recovers the usual plane-wave solutions ipR t L 
The BC (fTBl) imply: 


A R,Le 


.^fiEx 


Ar exp 
Ar exp 

which, eliminating Arl, yields 


lE 




= —i Al exp 
= i Al exp 




%E 




exp 


ME 


(e hL - l) 


= - 1 . 


The eigenmodes are then given by 

j-, _ h,Tr(m + 1/2) 
m “ 2(e hL - 1) 


. 7r (m + 1/2) 

= a{z) y 2L -, m = 0,±1,... 


( 21 ) 


( 22 ) 


(23) 


where a (z), defined as 


cl{z) = 


e z - 1 


^ = hL. 


(24) 


was interpreted in reference |ll] as the Fermi velocity. Notice that z has a finite value in the 
continuum limit since it can be written as z = ( h/a)C. In the latter reference it was shown that 
the single-particle spectrum of the Hamiltonian ([7|) corresponds to Eq. (1231 ) with a function a(z) 
given numerically in Fig. [21 As we can see, the analytic expression (1241) [22] gives a very good fit 
of the numerical data. 



Figure 2. The scaling function a(z) that gives the spectrum of the Hamiltonian (0 near the Fermi energy. 
The points are the numerical data, i.e., the slope of the spectrum at the Fermi point, and the continuous 
line is the analytic result. 
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Figure 3. Wavefunctions just below the Fermi level (m = 0, left column) and the next below the Fermi level 
(to = 1, right column) for 2 = 1 (top) and 2 (bottom). The analytic values of are given by Eq. (EUl) . 

To find the eigenfunctions with energy E m , we first compute the constants Ar^l using (ED 

Ar.l = exp j±'i ^ (e hL - l) - j J , (25) 

and Eq. m, obtaining 


V4 m) ^e fc|n| / 2 cos 


7T (n — m) . , 7r(m + 1/2) — 1 

— 1 -- + sign(n)—--—— —tt- 

2 B K ’ 2 e hL - 1 


(26) 


where n = ■ ■ ■, ± (^L — ^ and m = 0, ±1, • ■ - . Fig. [3] shows the numerical and analytic values 

of f° r m = 0,1 and 2 = 1 and 2. As hn = zn/L, we see that for the same value of 2 , all the 
curves collapse when expressed in the scaled variable n/L. 

The results obtained so far suggest that the continuum Hamiltonian (11711 can be brought to the 
standard canonical form of a free fermion with OBC. To show that this is indeed the case, let us 
make the change of variables 

e h\x\ _ i 


x = sign(x)- 


that maps the interval x € [— L, L] into the interval x € [— L, L] where 


(27) 


L = 


e hL -l 


(28) 
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The fermion fields in the variable x are given by 


dx W 1 / 2 


^R, L (x) = 1pR,L(x) = e h ^ /2 ^R,L(x), 


(29) 


that plugged into (TT7T) gives (recall that we set a = 1, so C = V) 

rL 


H ~ i 


dx 


l-L 


^ R dx4’R - V L dx'4’L 


(30) 


That is just the free fermion Hamiltonian for a chain of length 2 L. This result suggests that one 
could try to derive some of the properties of the rainbow Hamiltonian, Eq. CD, from those of 
the free fermion system. This will be done in the next section when discussing the entanglement 
properties of the GS. 

Notice that Eq. (12711 is not analytic at x = 0, but if we take x > 0 we obtain 

e hx - 1 


x = 


h 


(31) 


which is a conformal transformation (similarly, x = —(e~ hx — l)/h if x < 0). If we add the euclidean 
time coordinate, that is, x —>• x + ir. the transformation (1311) becomes periodic in r with a period 
equal to (3 = 2n/h. This result leads us to associate to the system an effective temperature 


h 

T = —. 
2vr 


(32) 


This result will be interpreted below. 


A. Validity of the continuum approximation 

As can be seen in Figs. [2]and[3l the continuum approximation provides very accurate predictions 
regarding the wavefunctions near the Fermi point and the Fermi velocity. Indeed, that is the 
expected range of validity of any continuum limit: the long-distance physics which takes place near 
the Fermi point. 

We have explored the limits of the validity of the continuum approximation. Fig. [I] shows 
the overlap between the predicted and the numerical single-particle wavefunctions as we go deeper 
beneath the Fermi surface. The horizontal axis shows the rescaled wavefunction index m/L , which 
is 0 for the Fermi level and 1 for the deepest one. The vertical axis corresponds to the overlap 
between the continuum approximation and the actual wavefunction, defined as 


Q — (Ip cord | 1pexact)\ ■ 


(33) 
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Figure 4. Top left: overlap between continuum approximation and numerical single-particle wavefunctions, 
in a one-to-one basis, for 2 = 1. Top right: overlap between continuum and numerical many-body states, as 
a function of z, for different values of L. Bottom: validity region for the continuum approximation, showing 
the lines of 90% and 95% overlap between the continuum and the exact wavefunctions. 

The numerical experiments were performed for L in the range of 50 to 500 and z = 1. Notice 
that for small wavefunction index m/L, the overlap is virtually one, but it decreases very fast 
behind a certain critical value. The explanation is that single-body wavefunctions which are deep 
below the Fermi energy vary over very small length scales, rendering the continuum approximation 
inaccurate. 

Even if the wavefunctions are not correctly predicted in a one-to-one basis, the complete Slater 
determinant composing the states can be similar. This possibility is checked in the top-right panel 
of Fig. U where we plot the overlap between the full Slater determinant states (continuum limit 
and numerical computation) as a function of z for different values of the system size L. We can 
see that below a certain critical 2 , the overlap stays close to one, and then it decreases to zero. 

The bottom panel of Fig. [5] shows the region of validity of the continuum approximation in 
the (L, z ) plane, by depicting two lines which mark the level 0.95 and the level 0.9 for the overlap 
between the numerical GS of the Hamiltonian and the continuum approximation obtained by the 
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Figure 5. Comparing the half-chain entropy for different L and h with the theoretical prediction, Eq. (1351) . 


deformed uniform wavefunctions. 


IV. ENTANGLEMENT STRUCTURE 


The most relevant fact about the entanglement structure of the GS of Eq. (ED is that the von 
Neumann entropy scales linearly with the system size —i.e., with a volume law- as soon as a < 1, 
with a prefactor that was determined numerically in [ll| as ~ — log(a)/6. Can this prefactor be 
explained? 


A. Entanglement in the continuum approximation 

The von Neumann entropy of the left half of a critical system with open boundary conditions 
and size L is given by 


Scft(L) = - log(L) + c', 
b 


(34) 


where c is the central charge of the model, and c' is an additive constant that includes the boundary 
entropy plus non-universal contributions [6, 7]. For the free fermionic system under study, a = 1, 
we have c = 1. Taking in combination Eqs. (1341) and (1281) . we can provide a prediction for the 
entropy of the half-chain in the deformed GS of (|7|). Indeed, substituting L by L, we obtain 


c ( — 1 

Scsp(L) = -log —— 


+ c , 


(35) 


which is checked in Fig. [5] (left) for low values of h, although its validity ranges far beyond that 
regime close to the conformal point. 
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Expression (1351) can be expanded in the limit when hL is large enough as 


Scsp{L ) « -hL, 
b 


in agreement with numerical estimation (15) of 



(36) 


Sl ~ —0.318 L log(a), 


(37) 


since h/2 = — log (a). Notice also that in the limit h —> 0 + , one recovers the expression (1341) . 

It is worth to compare Eq. (1351) with the entropy of a thermal state at inverse temperature 
= 1/T in a CFT [7] 

Scft(L) « |log ((binh (y)) « (38) 

where we have taken the limit L 3> /3, which leads to an extensive entropy. Comparing Eq. (1381) 
and Eq. (1361) we obtain that 


T = 


1 

0 


A 

2tt’ 


(39) 


in agreement with Eq. (1321) . which was based on the analytic extension of the transformation 
employed to derive the continuum limit, Eq. (|27|) . In other words, we can assert that the rainbow 
state is similar to a thermal state with temperature given by h. 


B. Renyi entropies 


Let us focus our attention on Renyi entropies, defined as: 

S’A = Y“^ logTl7l ^4’ 


(40) 


where A is a block with £ sites and pa is the corresponding density matrix. The von Neumann 
entropy can be obtained as the limit n —» 1 of S^ \ The expression of S'A for the GS of the free 
fermion model on an open chain of length 2 L is given by 


5, 


(n) 


AA 106 


4 L 


■ sin 


7T 


7 t£ \ 

2 LJ_ 


+ c' n + fn COS (tt£) 


8 L fi t£\ 

— sm —■ 

7T \2 L ) 


-K/n 


(41) 


where £ is the length of the block A that is located at either boundary of the chain. The first term 
is the familiar CFT contribution with c = 1, while the second term contains the fluctuations at 
Fermi momentum kp = it/ 2. K is the Luttinger parameter which in our case is equal to 1. The 
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Figure 6. Renyi entanglement entropy at half chain. All panels show z = {0, 0.4, 0.8, • • •, 4} from bottom 
to top, and fits are made to expression (l44)l . From top to bottom, and left to right, orders n = 1 (von 
Neumann), 2, 3 and 4 of the Renyi entropy. 


constants f n 


have also been computed analytically in reference 



/i = -l 


fn = 


(1 + -L) 

\ 2 T 2n) 


(n > 1). 


(42) 


1 - n r (i + 

1 \2 ^ 2nJ 

Let us consider the left-half block. According to (1411) . the Renyi entropy of order n is given by 




12 


n 


4 L\ 


— 1 + - log - ) +c' n + fn C0S(7 tL) - 


7r ) 


8 L\~ x l n 


7r J 


(43) 


Moving into the rainbow phase, we can give a first estimate of the Renyi entropies using the 
SDRG prescription. According to it, they are all equal among themselves, and equal to the von 
Neumann entropy Jr]]. This approximation becomes exact only in the a —>• 0 + limit. Otherwise, 


we should make use of the following exact diagonalization strategy 


fl 24]: 


1. Obtain the occupied single-body wavefunctions if k . 

2. Compute the correlation matrix, C{j = (GS\c\cj\GS) = with i and j both inside 


the considered block. 
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3. Diagonalize matrix C and obtain its eigenvalues {v p }. 

4. The Renyi entropy is then given by S [g = y4-_ J2 P log (v™ + (1 — v p ) n )- 


The numerical computations performed as above can be compared to a natural extension of 
expression 63D: 

Cn{z) 


S l\ Z ) = (l + -) log ( —) + d n (z) + f n (z)(-l) L ( 

□ 


8 Ly 1,n 


(44) 


This equation is a generalization of the Ansatz made in reference LUI for the von Neumann entropy 
of the half-chain S^\z). The comparison is performed in Fig. [HI which shows the Renyi entropies 
at half chain for different values of z in each panel, fitting the parameters c n (z), d n (z ) and f n {z)- 
The Luttinger constant is kept as K = 1. Oscillations in all cases decrease as z increases, but they 
always increase with the Renyi order n. 

The functions c n (z), d n (z ) and f n (z), are shown in Fig. [Tj Their expression can be derived 
replacing L by L in Eq. ([2811 . and writing 


L — —h. 


(45) 


which yields 

Cn(z) = 1, d n (z) = c ' n + ^ (l + ^ log , fn(z) = fn (^—^j , (46) 

Fig. [7] shows the fitting coefficients for different orders of the Renyi entropy, Eq. (1441) . for 
systems of size 10 2 < L < 10 3 and for a range of values 0 < z < 20. Panel (a) shows the small 
variation (< 4 • 10^ 2 ) for c n (z) in all range of z. Panels (b-c) show d n (z) and f n {z), solid lines 
are given by Eq. (06|). Notice the perfect agreement between these expressions and the numerical 
results. 


C. Entanglement spectrum and thermofield states 

The entanglement spectrum (EE) of this model was already considered in |11]. It is given by the 
eigenvalues of the entanglement Hamiltonian He , that is defined in terms of the reduced density 
matrix pa as pa = e~ Hs . Let us summarize the results obtained in [11] and provide an interpre¬ 
tation from the perspective obtained in the previous sections. The entanglement Hamiltonian He 
has the form 


He = Y1 e P b p b P + /o> 


(47) 
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Figure 7. Fitting parameters in Eq. (Hi)) as a function of z: (a) c n (z), where the lines are for visual guide 
only; (b) d n (z), where the continuous line is given in Eq. (H(T1) (c) f n (z), where the continuous line is given 
in Eq. (l46l) . 


where b p , are destruction and annihilation fermion operators, and e p are related to the eigenvalues 


up, u p 

of the correlation matrix in the block A (for details see 
are given approximately by 


Q). 


For large values of L, the energies e p 


e„ ~ A L p, p = 


±1 ±2 __ 
2 ’ 2 ’ ' ' ' ’ 2 


±^- L : even, 


0,±l,±2,...,±^fi L : odd. 


(48) 


The level spacing A l in turn was found in [ll] to be related to the half-chain entropy Sl as 

J2 


Sl 


7r 


3A l 


which, using Eq. (1281) and (1321) . implies 


. 27T 2 7 t{3 

Al ~Tl=T- 


(49) 


(50) 


Hence, the density matrix can be expressed as 


PA^e-Pnc™, H CF T = jJ2P b l b P > ( 51 ) 

P 

where Hcft is the CFT Hamiltonian of half of the chain. Thus, the single-body entanglement 
energies should fulfill, for different values of L and z, the following law: 
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-10 -5 0 5 10 

Eigenvalue Index 

Figure 8. Entanglement single-body energies e p for different values of L (60 to 160) and z (5 to L in each 
case), multiplied by z/(2tt 2 ). As predicted, they collapse to the diagonal line, following Eq. (15!Z]l . e p ~ 2?rp/z. 


which we can see confirmed in the results of Fig. [HJ 

We then arrive at the conclusion that the CSP state can be written as 


I i>csp) = X! e PEn/ 2 \n) L | n) R , 


(53) 


where |and \ti)l are orthonormal basis for the right and left pieces of the chain whose Hamil¬ 
tonians are isomorphic to He ft in Eq. m- A pure state of the form (1531) is called a thermofield 
state and has been employed in connection with black holes and the EPR=ER conjecture [25, 26], 
In these studies, T = 1//3 is the temperature of the black hole that can be expressed as a/(27r) 
where a is the acceleration of a Rindler observer. Looking at Eq. (13911 . we see that the constant h 
plays that role in our model. 

Numerical evidence for other cases where p ~ & - h cft was explored before using different 


scalings 


27]. 


V. TWO-DIMENSIONAL EXTENSION 

A natural question is: can the ID results be extended to 2D? In other terms, can we find a 
local 2D Hamiltonian whose GS violates maximally the area law? We shall next show that this is 
indeed possible in a rather simple way. 

Let us consider a 2L x 2L square lattice whose sites are labeled by X = (x, y) with x,y E 
{±1/2, ±3/2, • • • , ±(L/2 — 1/2)}. We define a hopping Hamiltonian of the form: 

H = - ^2 Jx,X'J x c X ' ± h.c. (54) 

(X,X') 
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Figure 9. Left: a small region of the 2 L x 2 L square lattice used to generate the 2D rainbow state. The 
nodes represent the sites, and the number attached to the link is the associated hopping amplitude, given by 
a' x \, with x the horizontal coordinate of the middle point. Right: 3D view of the structure of the hopping 
distribution. 


Table I. Fitting parameters for the 2D entropy function in equation (1551) 


a 

A(a) 

B(a) 

C(a) 

x 2 

i 

0 

0.234 

0.29 

3-10" 5 

0.95 0.0053 

0.0940 

0.77 

10" 9 

0.9 

0.0116 

0.0330 

0.87 

10 -9 

0.75 0.0307 

-0.0225 

0.85 

10" 9 

0.5 

0.0594 

-0.015 

0.70 

10 -9 


with Jx.x' = F ((X -\-X')/2) is only determined by the center of the segment joining points X 
and X'. In our case, we choose F(x,y) = cr x \ to resemble the ID analogue. Fig. [9] represents 
graphically a small region near the center. 

Fig. [TD] shows the entropy per unit length of a block composed of the left half of the system 
sl{q) = Sl(o)/L for different values of the coupling parameter a. The solid lines represent fits to 
an expression of the form 


sl{o) ~ A(a)L + B(a ) log(L) + C(a ), 


(55) 


where a non-zero linear term denotes a volumetric behavior of the entanglement 
logarithmic term is added in order to predict the correct behavior for a —>• 1“ 


entropy, and the 
28]. The fits can 

be seen in table |U Notice the low values of y 2 for a < 1 and the increase of the volume coefficient 
with a. 
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Figure 10. Entanglement entropy of the left half of the system, per unit length of the boundary, for different 
values of a for a 2D concentric singlet phase. 

VI. CONCLUSIONS AND FURTHER WORK 


In this work we have extended the previous studies of the concentric singlet state, or rainbow 
state jiol . 11]. which is a deformation of a ID system which exhibits a volume law for the entangle¬ 
ment entropy. We focus on a free fermion realization in order to benefit from the exact solubility. 
The extension presented in this article is based on the application of field-theory methods. We 
show that, in the vicinity of the conformal model, the ground state can be described by a map 
that is the union of the conformal maps associated to each half of the chain. The corresponding 
conformal transformations further suggest the definition of a temperature that is proportional to 
the parameter controlling the decay of the hopping parameters. We show how this deformation 
accounts for the change in the dispersion relation, the single-particle wavefunctions in the vicinity 
of the Fermi point, and the half-chain von Neumann and Renyi entropies. The appearance of a 
volume law entropy is linked to the existence of an effective temperature for the GS that is finally 
identified with a thermofield state. This striking result points towards an unexpected connection 
with the theory of black holes and the emergence of space-time from entanglement. 

Finally, we show how to extend the rainbow Hamiltonian to several dimensions in a natural 
way, and check that the entanglement entropy of the two-dimensional analogue grows as the area 
of the block. 
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Figure 11. Qubistic picture of the rainbow state for different sizes of the system 2 L = 10 and 2 L = 12 (upper 
and lower rows) and three values of a = 0.01, 0.3 and 1 (left to right columns). Color intensity denotes the 
wavefunction amplitude and hue denotes phase: red is positive and green negative. 
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Appendix A: Qubistic picture 

Qubism {lb| is a pictographic representation for quantum many-body states with the peculiarity 
that it allows for the visualization of entanglement. In summary, an N qubit wavefunction is shown 
on the [0, l] 2 square divided into 2 N / 2 x 2 N / 2 cells. Each of the 2^ wavefunction component are 
depicted into one of the cells, following a recursive pattern, in which the i -th qubit is associated 
with the i-tli length scale, in decreasing order. 

Fig. |TT] represents the qubistic plots of the rainbow ground state for two different sizes, 2 L = 10 
and 2 L = 12, and three values of a = 0.01, 0.3 and 1. Therefore, the rightmost panels correspond 
to the ground state of the free fermion model, and the leftmost panels represent the rainbow states. 
Notice that the representation is formed only by a finite and small set of points. 

Entanglement between the first pair of qubits and the rest can be visualized in the following way 
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[16| . Break the full square into 2x2 square of half-size. Count the number of different (strictly, 
linearly independent) images among the small squares. That number is an upper bound for the 
Schmidt rank, which is a measure of entanglement. The same procedure can continue, for the block 
composed of the first four qubits, if we decompose the original square into a 4 x 4 grid. In our case, 
notice that the dots in each of the small squares form a similar but different pattern. In fact, the 
number of different (independent) images coincides with the number of squares, 4 for the first two 
qubits, 16 for the first four, etc. This shows that the Schmidt rank grows as 2^, i.e., entanglement 
is maximal. 
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